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The angular distribution of the intensity in photoemission experiments is affected by electron 
diffraction patterns and by a smoothly varying contribution originated by both intrumental details 
and physical properties of the samples. The peculiar variety of scattering configurations available at 
the ALOISA beamline experimental station in Trieste stimulated the developement of an analytical 
description for the smooth angular dependence sustaining the diffraction features. We present here 
the basic formulae and their application to experimental data taken on the Fe/Cu3Au(001) system 
in order to highlight the role of the various parameters included in the distribution function. A 
specific model for the surface illumination has been developed as well as the overlayer thickness and 
surface roughness have been considered. 
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I. INTRODUCTION 



The study of the angular distribution of the Auger- and photo-electrons at medium kinetic energies (a few hundreds 
of eV) is a well established method to determine the atomic structure in a local environment surrounding the emitting 
atoms close to the surface of ordered materials. The lattice geometry is usually retrieved exploiting the observation 
of intensity maxima along the direction of close packed atom rows. By taking into account the electron diffraction 
(ED) features, one can also reconstruct the atomic distances between the emitter and its closest neighbors. 

In general, the distribution of the intensity in the ED patterns is originated by two contributions: one anisotropic 
component x, which is determined by the geometry of the local atomic structure, and one, slowly varying, isotropic 
component ISO which depends from both instrumental factors (such as sample illumination and detector angular 
resolution) and material dependent factors (such as atomic differential cross section, film thickness/escape depth, 
surface morphology/roughness). 

The origin of the various contributions to the ISO component has been identified since many years . Some 
aspects relevant to photoelectron diffraction have been recently reviewed [Q , in particular regarding the role of dipolar 
emission in surface science experiments. Nonetheless in this work we present original developemcnt of the ED analysis, 
which arises from the evolution of instrumental performance, in terms of analyzers positioning and angular resolution, 
as well as collimation and size of X-ray beams in third generation synchrotron sources. The analytical treatement 
of the instrumental factors is presented in detail for the end station of the ALOISA beamline (Trieste Synchrotron) , 
where a wide variety of scattering geometries is available for ED experiments. 

In the following, we present a functional form of the ISO component using physically meaningful parameters, that 
can be fitted at once with the x component to rigorously deal with the ED analysis. This functional form has been 
coupled to multiple scattering calculations of the anisotropic structure-dependent x component (we made use of the 
MSCD code by Chen and Van Hove ||]). This fitting procedure was applied to the structural study of thin Fe films 
grown on the Cu3Au(001) surface, where the Fe structure changes from an fcc-like configuration to a bcc-like one as a 
function of the film thickness. ||6|7| For this system, Auger- and photo-electron ED data where taken in combination 
with in-plane X-ray diffraction. 8| The latter technique yielded the lateral lattice parameter of the growing film with 
utmost precision and its value was used as input in the fitting of the ED patterns, the vertical lattice spacing thus 
being the only structural parameter to be determined. In the next section, we have compared our model of the ISO 
function with the ED data taken for a few Fe films to enlight the weight of a proper choice of the ISO functional 
form in the determination of the x component. 

II. EXPERIMENTAL SETUP 

The experimental chamber of the ALOISA beamline hosts several electron and photon detectors in UHV for the 
study of surface photoemission. X-ray and photolectron diffraction. This experimental setup exploits the very wide 
photon energy range given by the ALOISA grating-crystal monochromator, which spans from the dominion of X-ray 
photoemission spectroscopy (hundreds of eV) to that of X-ray diffraction (thousands of eV) Q . The three angular 
degrees of freedom given by the manipulator (a modified CTPO manipulator by Vacuum Generators) allow one to 
freely orient the sample surface with respect to the photon beam direction and the beam linear polarization. At the 
same time, the combined rotation of the frame (rotation B) hosting in UHV the electron analyzers, and of the whole 
experimental chamber (rotation C) allows one to explore a wide portion of the hemisphere above the surface for any 
sample orientation, (see Fig. |l] for a sketch of the scattering geometries accessible by the detectors). The electron 
analyzers are of the hemispherical type with a 33 mm mean radius and an optical lens system yielding an acceptance 
angle of '--^ 1° (FWHM) and a field of view of ~ 1 x 4 mm^ [|o|. The transverse size (FWHM) of the photon beam at 
the sample position is ~ 20 /im vertically and 150 /im horizontally, with a slight dependence on the photon energy. 
The photon beam is linearly polarized, with the electric field vector E in the horizontal plane [Tl|] . 

ED polar scans are usually collected by rotating one electron analyzer in the scattering plane from the surface normal 
to the horizon (polar angle 9), even though different geometries can be suitable to specific experimental requirements. 



such as photoelectron holography |12| or near node photoelectron diffraction ||13|. As far as the data reported in this 
paper are concerned, the scans have been taken in transverse-magnetic polarization with the substrate [001] direction 
oriented along the photon beam direction for a few selected grazing angles a on the surface. In the present case, the 
ED patterns have been collected for the Fe L23M23M45 Auger line at a kinetic energy of 698 eV (we used a photon 
energy of 900 eV). In one case, we measured the ED pattern from a policrystalline Fe sample for the 3p core level at 
a kinetic energy of 803 eV. The signal was always taken at the maximum of the corresponding spectral line and at 
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two suitably chosen energies aside the peak, in order to allow an effective subtraction of the background of secondary 
electrons. The films were prepared by in situ evaporation of Fe on a Cu3Au(001) surface and the thickness was 
measured in situ by X-ray reflectivity (further detail about the sample preparation can be found elsewhere [||). After 
an initial pseudomorphic stage, the Fe films assume a tetragonally distorted bcc structure, displaying a (001) surface 
orientation, but rotated by 45° with respect to the substrate. At the largest thickness considered here (36 A), the Fe 
film has almost recovered its bcc bulk structure, with a lateral lattice spacing a — 2.830 ± 0.005 A, as measured by 
in-plane X-ray diffraction 



III. DESCRIPTION OF THE 750 FUNCTION 



a. Angles definition A schematic picture of the angular degrees of freedom available in the determination of the 
scattering geometry is shown in Fig.|l|. The manipulator rotation is used to select the azimuthal orientation of the 
surface, without altering the grazing angle which is selected by the rotation a. The surface normal can be rotated 
around the beam axis by ±180° with the manipulator rotation ^: at ^ = 0° the surface normal versor n lies in the 
plane normal to the electric field vector E. 

The rotation B of the electron analyzer frame sets the deflection angle between the X-ray beam and the photoelectron 
wavevector fc, with a clearence of ±115°. Furthermore, the frame rotation axis B can be oriented in the y — z plane of 
the laboratory reference system exploiting the C rotation of the whole experimental chamber (independently of the 
manipulator settings). In particular, B is related to the rotation C as B — zcosC — ysmC. 

According to these definitions, the polar emission angle 6 (defined by k and the surface normal direction n) is 
defined as: 

cos 9{£^,a, B,C) — — cos B sin a + sin B cos C cos a sin ^ + sin B sin C cos a cos ^, (1) 

which simplifies to 6 = 90° — (B — a) if the analyzer is rotated in the plane containing the surface normal, which 
is often the case if the projection of the photoelectron wavevector onto the surface plane is to be mantained along a 
specific azimuthal direction, as in ED polar scans. 

In the following we describe the factors contributing to the measured ISO electron yield. 

b. Atomic cross-section In the dipole approximation (we will not consider multipolar expansion coefficients in 
the following applications), the atomic differential cross-section dani/d^ for the nl initial state of the photoemission 
process takes the well known analytical expression, 

^(/3;7)«[1 + 1/3(3 cos^-1)], (2) 

where /3(nZ, hv) is the asimmetry parameter in the matrix elements for the photoemission process jlj] and 7 is the 
angle between the directions of the polarization vector E and the photoelectron wavevector k (see Fig.|l|) . The value 
of cos 7 = sinScosC is determined only by the B and C rotations of the electron analyzer, independently from the 
sample surface orientation. 

c. Escape depth The effect of inelastic scattering on the probability of escape of photoelectrons from the surface 
of condensed matter samples is described by means of the inelastic mean free path (IMFP), whose variation as a 
function of the electron kinetic energy is well known and can be evaluated by analytical formulas psf . 

The flux of electrons emitted at a depth z and detected at a polar angle 9 from the surface normal will be reduced 
by inelastic scattering according to the Beer-Lambert relationship I{z,9) = e~^^^ . 

Integration over the depth from the surface yields the angular dependence expected for an emitting slab of matter 
of thickness D 

IiMFpiD,X;9) ^ I{z,9)dz^\cos9[l-e-^] (3) 
Jo 

which corresponds to the well known cos{9) behaviour in the limit 13 —> 00 of an homogeneous semi-infinite emitting 
volume. 

The same description can be used to describe the reduction of the photoemission intensity caused by a non-emitting 
layer, i.e. with different chemical composition, above the emitting region of the sample. Considering an emitting layer 
of thickness D below a non-emitting one of thickness D' , the full formula can be written as: 

IiMFp{D, A, D', A'; 9)^X cos9 [1 - e~^^] v€i? (4) 
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Further, the limited penetration of the X-ray beam can be taken into account: this yields a very small correction, 
since the IMFP of the photoelectrons at the kinetic energies typical of photoemission experiments (50-1000 eV) is of 
the order of 10 A, while the characteristic penetration depth of the X-ray photons is normally larger than 50 A. 

d. Surface roughness The angular distribution of the photoelectron intensity, specially for high values of polar 
angle, is also affected by the surface roughness of the emitting layer. Models can be set up, to take into account the 
effects of shadowing of the photoelectron flux at grazing emission, where the details of the resulting formulae depend 
on the details of the surface morphology |^ . 

A simple statistical model, which can be used in a phenomenological way, has been described in Ref . [|l6| . The rough 
surface is described in terms of an isotropic stationary random function z(x,y) of the in-plane surface coordinates 
with a mean value < z >= 0. The distribution of heights is assumed to be normal, so that the probability density for 

z is fiz) = —7= — . 

After a statistical analysis of the average shadowing by the surface protrusions, an analytical expression Iii(5,9) 
for the dependence on the polar angle is derived. Details of the formulae can be found in Ref. jl6|. Qualitatively 
speaking, assumes a constant unit value near normal emission, and drops to zero for 9 — > 90°, the steepness of 
the decay being determined by the amount of surface roughness. The shadowing can be completely neglected up to 
a maximum take-off value of the polar angle arctan((5), so that the phenomenological parameter 5 can be used as an 
effective marker of the surface roughness. 

e. Surface illumination Due to the high angular resolution of the ALOISA analyzers (acceptance angle of 1°, 
FWHM), the angular smearing of the ISO function originated by the angular acceptance can be neglected, although 
it must be taken into consideration in the simulation of the ED pattern x- ttiG other hand, the interplay between 
the size of the field of view and the illuminated portion area of the surface must be carefully taken into account. All 
the experiments at ALOISA are performed at grazing incidence, i.e. the value of the grazing angle a of the photon 
beam with respect to the sample surface is always in the [0 — 10°] range. Further, a is often set to values as low as 
a few degrees in order to exploit the increase of the signal from the surface region when total reflection condition are 
satisfied. Even if the transverse width Vbeam of the focussed photon beam (FWHM) is narrow, the beam footprint 
on the surface is elongated in the beam direction. The width of the illuminated surface area r(a) = Fbeam/sina can 
easily exceed the width of the surface area representing the projection of the analyzer slits (1 x 4mm^), specially when 
working in very grazing conditions (ai„ 1°) and in transverse magnetic polarization, ^ = 90°, where r^eam has its 
maximum value of ~ 150 /im. This polarization setting is indeed the most frequently used, since it takes the nodal 
planes of the atomic cross-section close to the surface horizon. 

The beam footprint on the surface can be considered as an unidimensional Gaussian intensity distribution 

exp |~ 2r^(Q;) } (s^^ -f^iS- H)' since the transverse distribution is always completely integrated by the analyzer slits. 

The integration of the photoemission intensity over the portion of illuminated area of the sample surface then results 
in the intensity factor 



f ' f 

IlLL{Tbeam,a;B,C) = / cxp < 

J-L/2 I 



L/2 r 



2(r(a)/2.35)2 



dx (5) 



The width L of the portion of illuminated area, which overlaps to the projected area of the analyzer slits onto the 
sample surface, can be derived analytically, for any orientation of the sample surface and detector position. For a 
generic scattering configuration, L is a non-trivial function of the experimental angles B, C, ^ and a, and can be also 
limited by the finite sample size. The complete description of the general case is given in the Appendix. For simple 
polar scans, in which the analyzer is rotated in the scattering plane defined by the photon wavevector and the surface 
normal , L takes the minimum value between the sample size and the projected slit width on the surface A/ cos 9. 
(see Fig. |) 

/. Complete ISO formula The complete formula combines all the described factors: 

ISO^A- ^iP) ■ IiMFp{D/X, D'/X') ■ In{6) ■ //LL(r,,,„, a), (6) 
dll 

where, for the sake of clarity, only the fitting parameters have been explicitely reported in the arguments of each 
correction factor. The cross-section asymmetry parameter /3 and the IMFP A, A' are set to calculated values. A is a 
scale factor to be determined by fitting the data as well as the roughness parameter S. The X-ray beam width Theam 
and the grazing angle a are usually set to the nominal values in the first step of the iterative fitting procedure. The 
emitting layer thickness D and the non-emitting overlayer thickness D' are also used as fitting parameters, unless an 
independent measurement is available for them (for instance by means of X-ray specular reflectivity Q). 
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IV. DATA ANALYSIS 



In Fig. H a series of polar scans of the Fe Auger LMM line taken for a few Fe films on Cu3Au(001) are shown. For 
each film thickness, the experimental curve is compared with the calculated ISO distribution. 

In order to emphasize the dependence on the emitting thickness D (see Eq. the curves were divided by a 
cosO factor. Besides the structural transition from fcc(lOO) to bcc(100)i?45° (witnessed by the angular shift of the 
main forward focusing directions), the weight of the thickness of the emitting layer on the ISO shape is remarkable, 
particularly at the highest values of polar emission. The calculations were performed using the nominal Fe overlayer 
thicknesses (as determined by X-ray reflectivity) and an IMFP of 13.4 A. The roughness parameter 6 was used as a free 
parameter and adjusted with a minimization procedure to the experimental data. For all the considered thicknesses, 
it fell in the range 5 = 1.0 ± 0.1 (maximum un-shadowed takeoff angle ~ 45°). 

To enlight the role of the illumination correction factor we show here data taken from the surface of a polycrystalline 
Fe sample, where the anisotropy x component is absent and the ED pattern can be simply reproduced by the ISO 
component. The polar scans of the Fe 3p photoelectron intensity (taken at a kinetic energy of 803 eV) are shown in 
Fig. ^ in comparison to a few calculated ISO functions, built by successive inclusion of inelastic scattering, atomic 
cross-section and illumination correction factors. The complete ISO component accurately fits the experimental data 
with the same set of parameters for both the scans. The illumination factor I ill was calculated using Tteam = 150 /im 
and a set to the nominal values, as indicated in the Figure. The asymmetry parameter in the cross-section factor was 
set to (3 = 1.485, as calculated using one of the standard formulae |l^. As expected, the weight of the illumination 
factor on the ISO function is more significant at lower grazing angles. Small values of the grazing angle a are usually 
chosen in the experiment settings in order to exploit the enhancement of the surface signal, obtained when a is as 
small as the critical value for total external reflection (of the order of a few degrees for soft X-rays) . 

As far as the structural analysis is concerned, the x component is often extracted by substracting an ISO component, 
as obtained by interpolation of the ED pattern with a polynomial or a cosine function. The comparison between the 
so obtained experimental anysotropy and simulations based on structural models is then approached at a later stage. 

Here, we propose the direct fitting of the complete ED experimental pattern, where the simulated ED curve is 
built up using both the multiple scattering ED calculation and the ISO function, so that EDsim = ISO{\ -t- Xsim)- 
In Fig. I we present an example of full fit to a polar scan taken on the Fe Auger LMM line, for a 36A Fe film on 
Cu3Au(001). The data are presented as a function of the analyzer B angle (the normal emission direction, i.e. the 
zero for the polar angle 9, is indicated by the vertical line). The Xsim component was calculated for a structural 
model by assuming the lateral lattice constant a = 2.830 A, as measured by in-plane X-ray diffraction and the 
vertical parameter c to be determined by the fitting procedure. The EDsim was fitted to the experimental EDg^p, 
using in this case the roughness parameter 6 and the X-ray beam width F as fitting parameters. The layer thickness 
D was set to the value provided by the X-ray reflectivity calibration. In fact, being D already as large as SXjmfp, 
it would have yielded a contribution very similar to the D /lambda oo limit, strongly reducing the reliability of its 
determination by the minimization of the ED fitting procedure. An additional scale factor was also added to Xsim to 
take into account the ratio between the x E^nd ISO amplitudes. 

This procedure was iterated over a series of ED calculations, so that for each structure-dependent Xsirm the 
corresponding best fit ISO parameters can be found, thus eliminating any arbitrary a priori choice of the ISO shape. 
The best fit parameters of the ISO component were found to be J = 0.9 ± 0.1 and F = 160 ± 20 fim. 

Three additional curves are added in the Figure, as obtained by calculating the ISO component including only one 
factor each time. As anticipated before, the contribution due to the electrons escape depth is very close to the limit 
IiMFP ~ cos 9. At low emission angles, the ISO shape is strongly affected by the illumination and roughness factors. 
The cross-section contribution is not shown, since the Auger emission was assumed to be isotropic. In fact, an exact 
calculation would have required a detailed study of all the decay channels partecipating to the Auger process 0, 
but the dependence on the angular momenta of the decay channels decreases as the kinetic energy increases, so that 
it becomes almost negligible at the present value of 698 eV p^ . 

At the bottom of Fig.^, Xexp — {EDexp — I SO) /ISO and Xsim — {EDsim — I SO) /I SO are compared to highhght 
the quality of the fit. With the given lateral lattice spacing a = 2.830 A, we found the ratio of the vertical to lateral 
spacing to be c/a = 1.03 ± 0.02. The Fe film displays a structure very close to its bcc bulk one, a = 2.86 A and 
c/a = 1, the latter being recovered at a much higher thickness iQ. In general, this approach to the calculation of the 
ISO component leads to a substantial improvement of the reliabilty of the structural model (with an indetermination 
not exceeding a few percents). 
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V. APPENDIX 



The accurate description of the portion of illuminated area, which is included in the projected area of the analyzer 
slits onto the sample surface, is derived here for a generic scattering setup {B, C, a). For the sake of clarity, we will 
assume ^ = 0°, since only the mutual orientation of ^ and C is relevant. 

First, the orientation of the slits in the laboratory reference system is defined by the vectors Bsiu and C'snt which 
connect the center of the slit with the center of the short and the long sides of the slits, respectively. The B and C 
labels recall that the wide aperture is along the direction scanned by the B angle, while the narrow one corresponds 
to a movement of the C rotation. 



B. 



AB 



lit 



k X Bn 



AB 



(sin B, — cos B cos C, — cos B sin C) 



^ AC - - AC 
Csut = —j- k X Bsiit = (0,sinC,cosC), 



(7) 

(8) 



where k = (cos B,sin B cos C, sin B sin C) is the photoelectron waveversor, Baxis = (0, — sin C, cos C) is the 
orientation of the B rotation axis and AS and AC are the full slit aperture in the two directions. 

The displacement of the four vertices of the slit from the center of the slit in the laboratory reference system is then 



K- = ±BsHt ±CsHt; J ^1,2,3 A- 



(9) 



The projection of the slit on the sample surface is in general a parallelogram, whose vertices are in the following 
positions in the laboratory reference system 



n ■ k 



k + V,., 



(10) 



where h is the surface normal waversor. The corresponding S" positions in the (2-Dimensional) surface reference 
system are 



S'j — {Sj ■ xtiit,Sj ■ yuit), 
xtiit = (cos Q;,0,sin a), 

ytilt = VLab- 



(11) 

(12) 
(13) 



The width L(B,C, a) of the portion of illuminated area falling inside the parallelogram can be computed as the 
minimum among the 



S" - S" 

qi _ or ■''f 



s'. 



;j = 1,2,3,4; 



(14) 



which are the crossing points of the parallelogram sides with the y = plane on the surface (we recall that the X-ray 
spot on the sample can be considered as monodimensional) . The integration lenght on the surface may be further 
limited by the surface physical length. 



[1] J. Cooper and R.N. Zare, J. Chem. Phys. 48 (1968) 942. 

[2] C.S. Fadley, R.J. Baird, W. Siekhaus, T. Novakov, S.A.L. Bergstrom, J. Elect. Spectrosc. Relat. Phenom. 4 (1974) 93-137. 
[3] C.S. Fadley, The Study of Surface Structures by Photoelectron Diffraction and Auger Electron Diffraction in Synchrotron 

Radiation research: Advances in Surface and Interface Science, Vol. 1: Techniques, edited by R. Z. Bachrach, Plenum 

Press, New York, (1992), and references therein. 

[4] D.P. Woodruff, J. Electron Spectrosc. Relat. Phenom. 100 (1999) 259-272. 

[5] MSCD is freely distributed by Y. Chen and M.A. Van Hove at: " http://electron.lbl.gov/mscdpack/mscdpack.html " ; also 

see Y. Chen et al, Phys. Rev. B 58 (1998) 13121. 



6 



[6] M.-T. Lin, J. Shen, W. Kuch, H. Jenniches, M. Klaua, CM. Schneider and J. Kirschner, Phys. Rev. B 55 (1997) 5886; 
Surf. Sci. 410 (1998) 298. 

B. Schirmer, B. Feldmann and M. Wuttig, Phys. Rev. B 58 (1998) 4984. 

F. Bruno, S. Ter reni, L. Floreano, D . Cvetko, P. Luches, L. Mattera, A. Morgante, R. Moroni, A. Verdini, and M. Canepa, 



to be published, cond-mat/0103458 



L. Floreano, G. Naletto, D. Cvetko, R. Gotter, M. Malvezzi, L. Marassi, A. Morgante, A. Santaniello, A. Verdini, F. 
Tommasini and G. Tondello, Rev. Sci. Instrum. 70 (1999) 3855. 

R. Gotter, A. Ruocco, A. Morgante, D. Cvetko, L. Floreano, F. Tommasini, and G. Stefani, Nucl. Instrum. Methods A 
467-468 (2001) 1468. 

Details of the experimental chamber can be found at the web page www.tasc.infm.it/tasc/lds/aIoisa/aloisa.htmI. 

J. Wider, F. Baumberger, M. Sambi, R. Gotter, A. Verdini, F. Bruno, D. Cvetko, A. Morgante, T. Greber and J. 

Osterwalder, Phys. Rev. Lett. 86 (2001) 2337. 

A. Verdini, M. Sambi, F. Bruno, D. Cvetko, M. Delia Negra, R. Gotter, L. Floreano, A. Morgante, G.A. Rizzi, G. Granozzi, 
Surf. Rev. Lett. 6 (1999) 1201. 

J. J. Yeh, Atomic Calculation of Photoionization Cross-Sections and Asymmetry Parameters, Gordon and Breach Science 
Publ. (1993). 

S. Tanuma, C.J. Powell, D.R. Penn, Surf. Interface Anal. 20 (1993) 77. 
A.V. Yakovenko, J. Electron Spectrosc. Relat. Phenom. 74 (1995) 237-240. 

T. Greber, J. Osterwalder, D. Naumovic, A. Stuck, S. Hiifner, and L. Schlapbach, Phys. rev. Lett. 69 (1992) 1947. 
Y.U. Ydzerda, Surf. Rev. Lett. 4 (1997) 161. 

R. Rocow, C. Carbone, Th. Dodt, F.P. Johnen and E. Kisker, Phys. Rev. B 41 (1990) 3426. 



7 




FIG. 1. Degrees of freedom in the choice of the scattering 
geometry for the ALOISA experimental station. The rotat- 
ing electron analyzer is represented by its entrance slits ES; 
the projection of the slits onto the sample surface is drawn 
as a white parallelogram, defining the portion of surface con- 
tributing to the electron yield. The X-ray beam footprint on 
the surface (illuminated area) is shown as a shadowed ellipsis. 
The manipulator rotations a, (p and £, are used to select in- 
dependently the grazing angle, the azimuthal orientation and 
the orientation of the surface with respect to the beam polar- 
ization vector E. The analyzer rotations B and C allows one 
to survey most of the sky over the sample surface. The polar 
angle 9 (emission angle referred to the surface normal n) and 
the 7 angle (between the photoelectron wavevector k and E) 
are also shown. 
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FIG. 2. Schematic model of the intensity distribution 
of the photon beam illuminated area on the sample surface, 
for polar scans taken rotating the analyzer in the scatter- 
ing plane. The projection of the analyzer slit on the surface 
L{9) — A/ cos 9 provides the integration limits in Eq. ^ 




20 40 
Polar angle (0) 

FIG. 3. Polar scans taken for the Fe Auger LMM line 
from a growing crystalline Fe overlayer on Cu3Au(001). The 
thickness dependence of the ISO distribution, also shown in 
the picture, has been highlighted by dividing both the data 
and the ISO curves by cos 9. The thicknesses have been in- 
dependently determined by x-ray reflectivity. 
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FIG. 4. Polar scans of Fe 3p intensity (at a kinetic energy 
of 803 eV) from the surface of a policrystalline Fe sample 
are compared with calculated ISO functions, built by succes- 
sive inclusion of the different factors (with /3= 1.485 for the 
atomic cross section factor and Fteam^lSO fj,m for the illu- 
mination factor). The full ISO function fits the experimental 
data with the same set of parameters (see text) for both the 
scans, taken at different value of grazing angle a. The verti- 
cal full line indicates the surface normal direction, while the 
vertical dot-dashed line is the direction of the photon beam 
electric field. 
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FIG. 5. Fe Auger LMM polar scan taken for the 36 A Fe 
film on Cu3Au(001). The experimental ED pattern is shown 
together with the full simulation EDsim and the correponding 
ISO function. At the bottom, x^xp = {ED^^p - ISO) /ISO 
and Xsim ~ {EDsim — I SO) /ISO are compared. Three addi- 
tional curves are added, as obtained recalculating the ISO in- 
cluding only one factor at a time. The cross-section contribu- 
tion is not shown, since we approximated the Auger emission 
with an isotropic distribution. The normal emission direction 
is indicated by the vertical line. 
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